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Abstract 

We comment on the Lee- Yang zero analysis for the study of the phase struc- 
ture of QCD at high temperature and baryon number density by Monte-Carlo 
simulations. We find that the sign problem for non-zero density QCD induces 
a serious problem in the finite volume scaling analysis of the Lee- Yang zeros 
for the investigation of the order of the phase transition. If the sign problem 
occurs at large volume, the Lee- Yang zeros will always approach the real axis 
of the complex parameter plane in the thermodynamic limit. This implies 
that a scaling behavior which would suggest a crossover transition will not 
be obtained. To clarify this problem, we discuss the Lee- Yang zero analysis 
for SU(3) pure gauge theory as a simple example without the sign problem, 
and then consider the case of non-zero density QCD. It is suggested that the 
distribution of the Lee- Yang zeros in the complex parameter space obtained 
by each simulation could be more important information for the investigation 
of the critical endpoint in the (T, n q ) plane than the finite volume scaling 
behavior. 
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I. INTRODUCTION 



In the last few years, remarkable progress in exploring the QCD phase structure in the 
temperature (T) and quark chemical potential (/i q ) plane have been made in numerical 
studies of lattice QCD. The phase transition line, separating hadron phase and quark-gluon 
plasma phase, was investigated from \x q = to finite fi q [1-6], and the equation of state was 
also analyzed quantitatively at low density [4,7-9]. 

Among others, study of the endpoint of the first order phase transition line in the (T, /j, q ) 
plane, whose existence is suggested by phenomenological studies [10,11], is particularly im- 
portant both from the experimental and theoretical point of view. To locate the critical 
endpoint, Fodor and Katz [2,3] investigated the positions of the Lee- Yang zeros (to be ex- 
plained more below) in the complex (3 = Q/g 2 plane using lattices with different spatial 
volumes, and examined the finite-volume scaling behavior of a Lee- Yang zero closest to the 
real axis. There are also studies in which the behavior of the critical endpoint as a function of 
the quark masses is examined by using the property that a critical endpoint exists at /i q = 
in the very small quark mass region for QCD with three flavors having degenerate quark 
masses. For extrapolating the result to the case with physical quark masses, an approach on 
the basis of the Taylor expansion in terms of fJ. q /T [12] and that of the imaginary chemical 
potential [13,14] have been developed. Moreover, a study of phase-quenched finite density 
QCD, i.e. simulations with an isospin chemical potential, has been discussed in Ref. [15]. 
The radius of convergence in the framework of the Taylor expansion of the grand canonical 
potential can establish a lower bound of the location of the critical endpoint [7,9,16]. Also, 
the Glasgow method [17] is an interesting approach for the study of QCD at non-zero baryon 
density. 

In this paper, we focus our attention on the method of the Lee- Yang zero [18] applied 
to finite density QCD. The Lee- Yang zero analysis is a popular method that is used to 
investigate the order of phase transitions. In order to study the existence of singularities 
in the thermodynamic limit (infinite volume limit), Lee and Yang proposed the following 
approach: Because there are no thermodynamic singularities as long as the volume is finite, 
the partition function Z is always non-zero; it can develop zeros only in the infinite volume 
limit. However, if a real parameter of the model is extended into the complex parameter 
plane, a singularity, characterized by Z = 0, can appear in the complex parameter plane 
even in a finite volume. These zeros are called the Lee- Yang zeros. Therefore, one can find 
the position of a singularity by exploring the position of the Lee- Yang zero in the complex 
parameter plane as a function of volume, and extrapolating the position of the Lee- Yang zero 
in the thermodynamic limit. For a system with a first order phase transition, the position of 
the nearest Lee- Yang zero approaches to the real axis in inverse proportion to the volume. 
On the other hand, the Lee- Yang zeros do not reach the real axis for crossover transition, 
i.e. rapid change without any thermodynamic singularities. For QCD at non-zero baryon 
density, we expect a rapid crossover transition in the low density regime which changes into 
a first order phase transition beyond a critical value of the density. 

As we mentioned, Fodor and Katz have investigated the finite volume dependence of the 
Lee- Yang zeros in the complex f3 = Q/g 2 plane for various values of the chemical potential. 
To carry out such analysis, the reweighting method [19] is adopted, in which one performs 
simulations at \i q = 0, and then corrects for the modified Boltzmann weight in the measure- 
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ment of observables. In this case a famous problem arises for large [i q jT and large volume, 
which is called "sign problem". The sign problem is caused by complex phase fluctuations 
of the fermion determinant. In the region of small fJ, q /T, the phase fluctuations are not large 
and the sign problem is not serious. However, if the sign of the modification factor changes 
frequently during subsequent Monte-Carlo steps for large fx q /T, the statistical error becomes 
larger than the expectation values in general. 

We find that the sign problem induces a serious problem in the finite volume scaling 
analysis of the Lee- Yang zero, which is used by Fodor and Katz. For any non-zero \i q the 
normalized partition function calculated on the real axis with necessarily limited statistics 
will numerically always be consistent with zero once the volumes grow large. This means 
that the scaling behavior suggesting a crossover transition will not be obtained for the case 
with the sign problem, which is in contrast with the usual expectation. 

Before discussing the case of non-zero density QCD, we study in Sec. II, as a simple 
example, the Lee- Yang zeros in the complex ft plane for SU(3) pure gauge theory by analyzing 
data from Monte-Carlo simulations. This model has a first order phase transition [20] and 
simulations are much easier than for QCD at non-zero baryon density. Moreover, the pure 
gauge theory does, of course, not have a sign problem. Hence it is a good example to 
demonstrate how the Lee- Yang analysis works in the complex ft plane. In addition, it 
will become clear during this exercise that the complex phase fluctuations arising from the 
imaginary part of ft near Lee- Yang zeros are quite similar to those coming from the quark 
determinant where the sign problem exists for non-zero density QCD. The problem of the 
complex measure is reviewed in Sec. III. There we also comment on the reweighting method 
for the study at non-zero baryon density. In Sec. IV, we discuss a problem which arises 
when we apply the Lee- Yang zero analysis for non-zero baryon density QCD by using the 
reweighting technique, and consider possible other approaches in the framework of the Lee- 
Yang zero analysis for the investigation of the critical endpoint. Conclusions and discussions 
are given in Sec. V. 



In this section, we apply the method of Lee- Yang zeros to the SU(3) pure gauge theory 
(quenched QCD). 1 The phase transition of the SU(3) pure gauge theory is known to be of 
first order [20], which is expected from the corresponding Z(3) spin models. The pure gauge 
theory is controlled by only one parameter ft = 6/g 2 with the partition function, 



where P is an averaged plaquette P = (Y^ x ,n<u W / /!^ 1 ( x ))/(6^V S ite), and W^ 1 is 1 x 1 Wilson 
loop operator for the lattice size A^ sitc = N% x N T . Here, N c and N T are spatial and temporal 



A pioneering study has been done for lattices with N T = 2 in Ref. [21]. 



II. LEE- YANG ZERO FOR SU(3) PURE GAUGE THEORY 



A. General remarks 




(1) 
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extension of the lattice. We extend the real parameter (3 into the complex plane (Ate, Am), 
and determine the position of Lee- Yang zeros, at which 2(Ate, Am) = is satisfied, by 
numerical simulations. We use standard Monte-Carlo techniques; configurations {U^} are 
generated with the probability of the Boltzmann weight. The expectation value of an oper- 
ator 0[Ufj], (O), is then calculated by taking an average over the configurations. We expect 
that the position of the Lee- Yang zero (/3 Re , A° m ) approaches the real f3 axis in the infinite 
volume limit; with f$( m ~ 1/V = N~ 3 for a first order phase transition. 

In carrying out the above calculation two problems arise: One is that the Monte Carlo 
method is applicable only to the expectation values of physical quantities but not to the 
partition function itself. Another problem is that the measure is complex for a complex cou- 
pling /3, and hence we cannot apply the Monte Carlo method directly, since the probabilities 
(Boltzmann weights) must be real and positive. To avoid these problems, we introduce the 
normalized partition function Z norm together with the reweighting technique, 

J £>£/ e 6/3 Ro AT sito P 

= /g6i/3 Im AT site P\ 

\ / (/3r c ,0) 

Here AP = P — (P) and \exp(6if3i m N sitc (P))\ = 1. Because the denominator Z((3^ e ,0) 
is always finite for any finite volume, the position of J?(Ate,Am) = can be identified by 
analyzing Z noim . Although the partition function is not zero for Am = 0, it can be zero at 
some points in the (Ate, Am) plane, when the complex phase factor in Eq.(2) changes sign 
frequently on the generated configurations. For the determination of the critical point in 
the original theory, i.e. on the real (3 axis, the position of the nearest Lee- Yang zero should 
be investigated function of the volume V = N 3 . 

The mechanism that leads the occurrence of a Lee- Yang zero in Eq.(2) is quite similar to 
that which limits the applicability of the reweighting method for QCD with finite chemical 
potential [4,22]. We will discuss this in more detail in Sec. III. At a point for which 
the width of the probability distribution of 6/3i m N S i te P is smaller than 0(tt/2), the sign of 
the complex phase does not change. Therefore, the standard deviation of the plaquette 
distribution is required to be larger than 7r/(12/3 Im ./V site ) in the region where Lee- Yang zeros 
exist. Moreover, because the square of standard deviation is in proportion to the value 
of the plaquette susceptibility, the position of its maximum must agree with the position 
where Z norm becomes minimal as function of Ate for fixed Am- Hence, the real part, Ate, of 
the position of the nearest Lee- Yang zero must be consistent with the peak position of the 
plaquette susceptibility. The method to find a critical point from the location of Lee- Yang 
zeros thus is essentially the same as the method which determines a critical point through 
the location of the peak position of the susceptibility and its finite volume scaling. 

Here, it is instructive to introduce a probability distribution function for the plaquette, 
w(P), which is defined by 

W (P') = 1 J VU5(P' - p)e 6 P** N «* p , (3) 

where 5(x) is the delta function. Then, Eq.(2) can be rewritten as 



-^norm (Ate) Am) 



^(Ate,Am) 
^(Ate,0) 



6i/3i m V aite AP 



(/3r c ,0) 



(2) 



4 



J dPe m ^ N ^ Ap w(P) . (4) 



This means that the partition function Z norm as a function of 6fli m N site is obtained through 
a Fourier transformation of w(P). 

Using this equation, the relation between the scaling behavior of the Lee- Yang zeros in 
the infinite volume limit and the distribution function of the plaquette becomes clearer. In 
Monte-Carlo simulations, configurations with probabilities proportional to their Boltzmann 
weight are generated by a computer, and we obtain a distribution function of the plaquette 
from the histogram of the plaquette. The histogram has usually a Gaussian shape at a 
normal, non-critical point, but it deviates from the Gaussian form near a critical point, and 
attains a double peak shape at a first order transition point, corresponding to the coexistence 
of two phases. 

For the case of a non-singular point of /3r c or a crossover pseudo-critical point, where 
the distribution is expected to be a Gaussian function, the point of Z norm = does not exist 
except in the limit of (3i m N s n e — > oo or — oo, because the function which is obtained through 
a Fourier transformation of a Gaussian function again is a Gaussian function. Of course, 
results of numerical simulations have statistical errors, hence Z norm can become zero "within 
errors", if the expectation value and the error become of the same order. However, in this 
case, the point at which Z norm = appears at random in terms of [3i m x iV S i te . Therefore, 
the volume dependence of the position of the Lee- Yang zero {0^, Pi m ) does not necessary 
to be $ m ~ l/V(= N- 3 ) for fixed JV T . 

On the other hand, in the case of a first order phase transition, we expect that the 
plaquette histogram has two peaks having the same peak height at the transition point. 
Performing the Fourier transformation of such a double peaked function leads to a function 
which has zeros periodically. For example, a distribution function w(P) having two Gaussian 
peaks at AP = ±A leads to a normalized partition function Z norm which has zeros at 

This is mathematically the same as that for the interference experiment using a laser and 
a double-slit. The Lee- Yang zeros correspond to dark lines (destructive interference) and 
they appear periodically as given in Eq.(5). Moreover, for a first order phase transition the 
difference of plaquette values in cold and hot phases, 2A, is related to the latent heat As, 
i.e. the energy difference between the hot and cold phases, 

where a is the lattice spacing. Since Ae is non-zero, A does not vanish in the infinite volume 
limit (N% — > oo). Therefore, we find that in the infinite volume limit the nearest Lee- Yang 
zero approaches the real (5 axis like /5° m ~ 1/V", which is consistent with the general argument 
on the Lee- Yang zero for a first order phase transition. We also emphasize that the isolated 
Lee- Yang zeros appear periodically. The distances to these points from the real axis are 
1, 3, 5, • • • in units of the distance to the first Lee- Yang zero. This is also an important 
property, which is not observed for a crossover transition. 
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In addition, the discussion given for the plaquette distribution function can be extended 
to the analysis of fourth order Binder cumulants, 

B * = (7) 

which is an alternative to the method of Lee- Yang zeros often used to identify the order 
of a phase transition. The value of the Binder cumulant at the critical point depends on 
the universality class. In the case of a first order phase transition, assuming the plaquette 
distribution is a double peaked function, the Binder cumulants are estimated as 

JdPAP 4 w(P) A 4 

where the distance between two peaks is 2A and is wider than the width of each peak. On 
the other hand, when the distribution function can be modeled by a Gaussian function for 
a crossover transition, the Binder cumulants are given by 




Jx/n JdPAP 4 e- xAp2 , . , . 

(y/x/TT J dPAP 2 e-* AP2 ) 

In a region where a first order phase transition changes to a crossover, the Binder cumulant 
changes rapidly from one to three. We expect to find such a region for full QCD at high 
temperature and density. The value of the Binder cumulant at the endpoint of the first order 
transition line, which is of second order, is determined by the universality class. Hence, the 
plaquette distribution function plays an important role for both methods to identify the 
order of a phase transition. 



B. Numerical results 



We calculate the normalized partition function for SU(3) pure gauge theory to find Lee- 
Yang zeros in the complex (3 plane, using data for plaquettes obtained by QCDPAX in 
Ref. [23]. There are five data sets measured at the transition point for N T = 4 and 6. The 
spatial lattice sizes are 24 2 x 36 x 4, 12 2 x 24 x 4, 36 2 x 48 x 6, 24 3 x 6, and 20 3 x 6. O(10 6 ) 
configurations are available for the analysis of each data set. The reweighting technique is 
also used for the real (3 direction to analyze the Lee- Yang zeros in the complex (3 plane for a 
data set obtained at only one (3 point (Ale, Am) = (A, 0). The normalized partition function 
is given by 



-2norm(Ale) Am) — 



n Am AUe AP p 6/3 Re N Bite AP\ 



I (A),o) 



(10) 



Figure 1 shows the contour plot of Z noTni for the 24 2 x 36 x 4 lattice. The simulation point 
is A — 3.6492. In this definition, Z novm is normalized to be one on the real (3 axis. Circles 
at (/^ e ,A°m) ~ (5.6925,0.0021) and (5.6931,0.0056) are Lee- Yang zeros. Since the SU(3) 
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pure gauge theory has a first order phase transition, Lee- Yang zeros appear periodically. For 
this data, two clear peaks are visible in the plaquette histogram [23]. The distance between 
these two peaks is 2A « 0.003. The positions of the Lee- Yang zeros are consistent with 
Am ~ 7r/(12iV sitc A) « 0.002 and 37r/(12A sitc A) « 0.006, as given in Eq.(5). 

Above property is not seen so clearly for lattices having small N a and large N T . The 
position of the next-to-leading zero points of Z norm appear at random for the other data sets 
relative to the nearest zero point. The positions of Lee- Yang zeros are shown in Table I. We 
could not obtain clearly isolated Lee- Yang zeros for the lattices 24 3 x 6 and 20 3 x 6. The 
second nearest Lee- Yang zero to the real axis could be measured only for the 24 2 x 36 x 4 
lattice. The result on the 36 2 x 48 x 6 lattice (/3 = 5.8936) is also shown in Fig. 2. Only 
the nearest Lee- Yang zero is obtained clearly. The Lee- Yang zero becomes less clear as N T 
increases and N a decreases, hence simulations on lattices having large N a /N T seem to be 
necessary for the study of Lee- Yang zeros. 

The values for P° m V on 24 2 x 36 x 4 and 12 2 x 24 x 4 lattices are 43.9(5) and 42.0(6), 
respectively. These are roughly constant and suggest the scaling behavior of /3° m ~ 1/V for 
a first order phase transition. Also, in the previous study for lattices with N T = 2 [21], the 
1/V scaling behavior has been confirmed for N a = 6,8,10 and 12. However, for a more 
precise quantitative investigation that takes into account the errors, the spatial lattice size 
12 2 x 24 may not be large enough to check the 1/V scaling for N T = 4, since the difference 
of (3j m V is larger than the statistical error. We should fit the data obtained on more than 
two lattices by a curved function of 1/V, to confirm through a 1/V scaling analysis that the 
phase transition of the SU(3) pure gauge theory is first order. E.g. for the study of the SU(2) 
gauge-Higgs model [24], the following fitting functions have been used, Im.Ko(V) = k,q+CV~ u 
and lmK, (y) = Kq + CV^ 1 + DV~ 2 for a complex parameter k in the model with fitting 
parameters Kq, C, D, and v. 

III. COMPLEX PHASE FLUCTUATION AND SIGN PROBLEM 

As seen in the previous section, the investigation of Lee- Yang zeros in the complex j3 
plane seems to be useful to identify the order of phase transition. However, if we try to 
extend this analysis to full QCD at non-zero baryon density, a serious problem arises. This 
problem is closely related to the sign problem for finite density QCD, since the normalized 
partition function can be zero in the complex (5 plane due to fluctuations of the complex 
phase related to (3i m and also due to the complex phase from the quark determinant that 
causes the sign problem. Before discussing the Lee- Yang zero analysis for finite density 
QCD, we would like to review the sign problem briefly. 

The main difficulty for studies at finite baryon density is that the Boltzmann weight is 
complex if the chemical potential is non-zero. In this case the Monte-Carlo method is not 
applicable directly, since configurations cannot be generated with a complex probability. 
One approach to avoid this problem is the reweighting method. We perform simulations at 
[i = 0, and incorporate the remaining part of the correct Boltzmann weight for finite [i in 
the calculation of expectation values. Expectation values (O) at (j3,[i) are thus computed 
by a simulation at ((3, 0) using the following identity, 
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(o) iM = 



detM(/x)-lndet Af(0)) 



09,0) 



^ e JV f (lndetM(^)-ln detM(O))^ 



(11) 



where M is the quark matrix and Nf is the number of flavors (iVf/4 for staggered type quarks 
instead of Nf); // is a quark chemical potential in lattice units, i.e. = fi q a = /i q /(N T T), 
and n q is the quark chemical potential in physical units. This is the basic formula of the 
reweighting method. However, because lndetM(/x) is complex, the calculations of the nu- 
merator and denominator in Eq.(ll) becomes in practice increasingly more difficult for larger 
/i. We define the phase of the quark determinant 9 by (det M(fj 1 )) Nf ^ 4 = | det M(fi)\ Nf ^ 4: e ie 
for staggered type quarks. If the typical value of 9 becomes larger than ir/2, the real part of 
e %e (= cos#) changes its sign frequently. Eventually both the numerator and denominator of 
Eq.(ll) become smaller than their statistical errors and Eq.(ll) can no longer be evaluated. 
We call it the "sign problem". 

Here, the denominator of Eq.(ll), or simply (cos 9), is a good indicator for the occurrence 
of the sign problem. If this indicator is zero within statistical errors, Eq.(ll) cannot be 
computed. In the following we give an estimate for the value of the complex phase. Since 
the direct calculation of the quark determinant is difficult except for calculations on small 
lattices, we expand In det M(/i) in a Taylor series, 



lndet Af(//) 



lndetM(O) = 

n=l 



d n (In det M) 
d/i n 



nf 



(12) 



Then, we can easily separate it into real and imaginary parts because the even derivatives 
of In det M(/x) are real and the odd derivatives are purely imaginary [4]. The complex phase 
9 is given by 



n N { ™ T d 2n - l (\ndetM) a 2n ~ l 

e = TL Im " 



dji 



2n-l 



(2ra- 1)!' 



(13) 



for staggered type quarks at small /i. The Taylor expansion coefficients are rather easy to 
calculate by using the stochastic noise method. The comparison between the value of 9 with 
this approximation and the exact value has been done in Ref. [25], and the reliability and 
the application range have been discussed. 

We use data for the Taylor expansion coefficients obtained in Ref. [9]. The data were 
generated by using Symanzik-improved gauge and p4-improved staggered fermion actions. 
Coefficients up to 0(/x 5 ) have been calculated. Figure 3 shows the indicator (cos#) measured 
at /3 — {3,60,3.65, and 3.68}, for ma = 0.1, corresponding to T/T c = 0.90, 1.00, and 1.07, 
respectively, on a 16 3 x 4 lattice 2 . We also estimate the values of ji q /T = N T ji at which 



2 As mentioned in Ref. [4], in the calculation using the stochastic noise method, the error due to 
the finite number of noise vector (iVnoisc) is large for the calculation of (cos 9) with ./Vnoisc = 10. 
For the purpose of this study we increased the number of noise vector to N no \ se = 100. We checked 
that the difference between the results with N no \ sc = 50 and A^ no i se = 100 is about 10% for the 
calculation of the position at which (cos#) = 0.1. 
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(cos 9) = 0.1, 0.01 and 0.0. Results are given in Table II. The situation now is quite similar to 
the non-singular case of the normalized partition function in the complex (3 plane discussed 
in the previous section. This becomes even more apparent if we consider for simplicity 
only the first term in the expansion of 9 which is proportional to /i. Then (cos 9), /i and 
Im[9(mdet M)/dfj] correspond to the normalized partition function, Am and plaquette, 
respectively. 

Because the distribution of the complex phase 9 is almost of Gaussian shape, the indi- 
cator, (cos 9), decreases exponentially as \i increases, and it may cross zero at a point where 
the expectation value becomes smaller than the statistical error. Therefore, the points of 
(cos 9) = appear accidentally and the results given in Table II are unstable. Moreover, 
<9(ln det M) / d\x becomes larger as the volume increases, hence the indicator for the sign 
problem vanishes in the infinite volume limit for any non-zero //, which means that the 
range of applicability for the reweighting method approaches /i = in the infinite volume 
limit. 



IV. LEE- YANG ZERO ANALYSIS FOR FINITE DENSITY QCD 



The Lee- Yang zero analysis for finite density QCD has been performed by Fodor and 
Katz [2,3]. They measured the normalized partition function Z norm using the reweighting 
method, and determined the points where Z norm = as a function of spatial volume. The 
normalized partition function is defined by 



i (Ale) Am, A 1 ) = 



^(Ate, Am, A*) 



or 



i (Ale, Am, A 4 ) 



^(Ale,0,0) 



e 6i/?i m iV s iteAP e i 



-^(Ale, Am, A 1 ) 



3 ( N { /4) (In det M(/i) -In det M(0)) 



(fee, 0,0) 



(14) 



£(Ale,0,AO 

^iPlmNsiteAP 



(Nf /4) (In det M(ji) -In det M(0)) 



(Ale ,0,0) 



l^ie | e (JV f /4)(lndetM(/i)-lndet Af(0))| 



'(fee, 0,0) 



(15) 



for staggered type quarks. 9 is the complex phase of exp[(7Vf/4)(mdet M(/x) — In det M(0))]. 
Since the numerator of Eq.(15), which is the same as Z norm in Eq.(14), is required to be zero 
at a zero point of Z norm in Eq.(15), we consider Eq.(14) as an indicator for the Lee- Yang 
zero. 

Here, we notice that for Am — this normalized partition function is exactly the same 
as the indicator for the sign problem, i.e. the denominator of Eq.(ll). As discussed in 
the previous section, in any practical simulation this indicator will be consistent with zero 
within errors for large values of ji. Moreover, the region where the indicator is non-zero 
becomes narrower as the volume increases, and this region vanishes in the infinite volume 
limit. Hence, the Lee- Yang zeros always approach the Am = axis in the infinite volume 
limit for any finite [i. This means that the scaling behavior for crossover will not be obtained 
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for the case with the sign problem. This is clearly different from the usual expectation for 
the QCD phase diagram in the (T, //) plane. Most model calculations suggest that the 
transition is crossover in the low density region. This might be a problem of the definition 
of the normalized partition function, Eq.(14). The normalized partition function on the 
real (3 axis is exactly one for /i — 0, but it vanishes for finite \i in the infinite volume 
limit. Therefore, it is very difficult to distinguish the first order transition and crossover 
by investigating the position of Lee- Yang zeros as a function of spatial volume. This is the 
most important difference between the definitions for the pure gauge theory in Sec. II and 
QCD at non-zero density. 

The critical endpoint is shown to be located at fi B — 3fi q = 725(35)MeV in Ref. [2] 
and 360(40)MeV in Ref. [3], which is inconsistent with the above argument. This may be a 
problem of the fitting function. In Refs. [2,3] the position of Lee- Yang zeros has been fitted 
by Pim = ^(V^O + Am; "where A and /3f£ are fitting parameters. The first order transition 
and crossover have then been distinguished by the value of (3^. As we discussed in Sec. II, 
this fitting function is too simple to fit the data of Lee- Yang zeros obtained on lattices as 
small as those used in Ref. [2,3], i.e. V < 12 3 . In fact, if one assumes a curved extrapolation 
function, all data in Table 1 of Ref. [3] seems to approach f3i m = in the 1/V — > limit. 

In our argument, the statistical error of Z norm , which is controlled by the number of 
configurations in the Monte-Carlo simulation, plays an important role. If the statistical error 
of Znorm becomes much smaller than the mean value of Z norm by increasing the statistics for 
each simulation, the method in Ref. [2,3] would be applicable. However, one cannot satisfy 
this condition in general simply because we are looking for the Lee- Yang zero which gives 
Z norm = 0. Namely, statistical error cannot be smaller than 0. Moreover, if the error of 
Z-norm is sizeable, there appear fake Lee- Yang zeros which are located even closer to the real 
(3 axis than the true zero in the region where the mean value of Z norm is smaller than the 
error. Since we adopt the closest zero as the Lee- Yang zero in the actual scaling analysis, 
we may thus misidentify the true zero by the fake one. 

The above point can be seen explicitly for the second Lee- Yang zero of the SU(3) gauge 
theory with N T = 6 shown in Fig. 2. Theoretically, we expect the second Lee- Yang zero 
exists around f3\ m « 0.013, i.e. three times larger than that of the first Lee- Yang zero as 
shown in Fig. 1. However we find several f3 which give Z norm = in the region of fli m > 0.006, 
and they distribute randomly. If we identify the second nearest point as the second Lee- Yang 
zero, the resulting f3i m is much smaller than the theoretical expectation. This problem is 
caused by the existence of the region in the complex /3-plane where the statistical error of 
■^norm is larger than the mean value as discussed above. As the statistics is increased for 
fixed V, such a region should become smaller and fake Lee- Yang zeros should disappear. 

Now we discuss how the above situation changes by increasing the volume V. Fortunately, 
in the SU(3) gauge theory, the location of Lee- Yang zero can be determined better also as 
the volume increases as shown in Sec. II. In this case, the validity of the scaling assumption 
becomes better as the volume increases. On the other hand, for the case with the sign 
problem, the error normalized by the mean value grows exponentially as a function of volume. 
Then the size of the region having fake Lee- Yang zeros cannot be made smaller unless one 
has exponentially large statistics. This leads to the conclusion that the quality of the scaling 
analysis is not improved by increasing the volume, and any reliable information about Lee- 
Yang zeros in the infinite volume limit cannot be obtained. This is the reason why the 
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serious problem of identifying the critical endpoint is intimately related to the sign problem. 

This discussion suggests that the finite volume scaling analysis suffers serious damage 
through the unsolved sign problem, and it is very difficult to apply the criterion used by 
Fodor and Katz for the investigation of the critical endpoint in practice. However, the 
property of the second nearest Lee- Yang zero characteristic for a first order transition in 
Fig. 1, i.e. the fact that the distance to the second Lee- Yang zero from the real axis is 
three times larger than that of the first Lee- Yang zero, could be investigated on a finite 
lattice. This study is possible for small \x without taking the infinite volume limit. On the 
other hand, we do not expect any isolated Lee- Yang zero for a crossover transition, hence we 
may be able to determine the order of phase transition by investigating the distribution of 
Lee- Yang zeros in the complex (3 plane. Although the measurements of the second Lee- Yang 
zero may require large lattice sizes and high statistics, as seen for the case of pure SU(3) 
gauge theory, it may be possible to find the region of the first order phase transition, if the 
critical endpoint exists in the low density region. 

V. CONCLUSIONS 

We commented on the Lee- Yang zero analysis for the study of the critical endpoint in 
the (T, phase diagram. It is found that the Lee- Yang zero analysis at non-zero baryon 
density encounters a serious problem. The complex phases of the quark determinant and the 
complex f3 are mixed at non-zero chemical potential. In this case, in practical simulations 
with limited statistics the normalized partition function can develop zeros even on the real 
f3 axis for large fi q in finite volumes. Moreover, in the infinite volume the normalized 
partition function is always zero except for fi q = 0. This means that the nearest Lee- Yang 
zero always approaches the real /3 axis in the infinite volume limit. The scaling behavior 
suggesting a crossover transition thus will not be obtained. This is clearly different from 
usual expectations for the QCD phase diagram. To avoid this problem, the sign problem 
must be removed by careful treatments increasing the number of configurations exponentially 
as the volume or fi q /T increases, otherwise the finite volume scaling behavior for the position 
of Lee- Yang zeros, which has been analyzed by Fodor and Katz [2,3], does not provide an 
appropriate criterion for the investigation of the order of the phase transition. 

To make the underlying problem more transparent, we applied the Lee- Yang zero analysis 
to the SU(3) pure gauge theory, which does not have a sign problem and for which the 
simulations are much easier. Lee- Yang zeros are found in the complex f3 plane. They appear 
periodically as expected by the discussion using a plaquette distribution function for a first 
order phase transition. The positions of the first Lee- Yang zero on two lattices having 
different volume sizes are roughly consistent with the finite size scaling behavior for a first 
order phase transition, i.e. /3° m ~ 1/V. However, for quantitative analysis it is necessary to 
fit data from more than two different lattice sizes by a curved function to study the order 
of the phase transition. It is found, in this analysis, that complex phase fluctuations arising 
from the imaginary part of (3 play an important role, and the mechanism that leads to the 
appearance of the Lee- Yang zeros is quite similar to the situation in QCD where the sign 
problem is present. 

The property of a first order phase transition that isolated Lee- Yang zeros appear peri- 
odically at j3i m ~ C(2n + 1), where C is the distance to the nearest Lee- Yang zero and n 
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is an integer, is free from the problems that arise in the infinite volume limit. Therefore, 
to investigate the pattern of the appearance of Lee- Yang zeros in the (/9R e , Am) plane is im- 
portant. For this calculation, the simulations by using high statistics data and large lattice 
size are indispensable. Further studies are clearly important to find the endpoint of the first 
order phase transition line in the (T, /i q ) plane. 

Recently, a close relation between the strength of the sign problem and the position of 
the phase transition line for pion condensation in phase quenched QCD has been discussed 
in Ref. [26]. There it has been found that the endpoints of the first order transition line 
determined in Ref. [2,3] are located near the phase transition line of pion condensation. 
These results may relate to our discussion given here. 

Moreover, the pathologies in the Glasgow method have been discussed in Ref. [17]. Sim- 
ilar problems arise also in the Glasgow method. It would be interesting to consider the 
relation between the pathologies in the Glasgow method and those in the Lee- Yang zero 
analysis. 
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TABLES 



TABLE I. Positions of Lee- Yang zeros for the SU(3) pure gauge theory 



Lattice size 






Am 


12 2 x 24 x 4 


1st zero 


5.69178(23) 


0.01216(17) 


24 2 x 36 x 4 


1st zero 


5.69252(5) 


0.00212(2) 


24 2 x 36 x 4 


2nd zero 


5.69309(7) 


0.00556(7) 


36 2 x 48 x 6 


1st zero 


5.89411(10) 


0.00434(8) 



TABLE II. Values of n q /T = N T fi at which {cob 9) = 0.1, 0.01 and 0. 7V site = 16 3 x 4. 



T/Tq (cos 9) = 0.1 (cosfl) = 0.01 (cosfl) = 0.0 

~O90 0.70(2) 1.0(2) 1.2(2) 

0.96 0.80(2) 1.1(1) 1.1(1) 

1.00 0.87(2) 1.9(2) 2.3(4) 

1.02 0.96(3) 2.2(12) 2.3(1) 

1.07 1.13(3) 1.8(4) 2.0(2) 
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FIGURES 



24 x36x4 lattice 
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FIG. 1. Contour plot of the normalized partition function Z norm in the (/?R e , Am) plane mea- 
sured on the 24 2 x 36 x 4 lattice. Values in the right edge are Z norrn . The simulation point is 
fa = 5.6925. 
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FIG. 2. Contour plot of the normalized partition function Z Dorm for the 36 2 x 48 x 6 lattice. 
The simulation point is (3q = 5.8936. 
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